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When galaxy clusters collide, they generate shock fronts in the hot intracluster medium. 
Observations of these shocks can provide valuable information on the merger dynamics 
^ and physical conditions in the cluster plasma, and even help constrain the nature of 

Odark matter. To study shock fronts, one needs an X-ray telescope with high angular 
resolution (such as Chandra), and be lucky to see the merger from the right angle and 
qq at the right moment. As of this writing, only a handful of merger shock fronts have been 

I discovered and confirmed using both X-ray imaging and gas temperature data — those 

in IE 0657-56, A520, A754, and two fronts in A2146. A few more are probable shocks 
l | awaiting temperature profile confirmation — those in A521, RXJ 1314-25, A3667, A2744, 

and Coma. The highest Mach number is 3 in IE 0657-56, while the rest has M ~ 1.6 — 2. 

U Interestingly, all these relatively weak X-ray shocks coincide with sharp edges in their 
host cluster's synchrotron radio halos (except in A3667, where it coincides with the 
^ distinct radio relic, and A2146, which does not have radio data yet). This is contrary to 

f*) ^ the common wisdom that weak shocks are inefficient particle accelerators, and may shed 

I light on the mechanisms of relativistic electron production in astrophysical plasmas. 

O 

Keywords: Galaxy clusters; Intergalactic medium; Shock fronts; Dark matter; Cosmic 
1. Introduction 

Mergers of galaxy clusters provide a unique laboratory to study the cluster physics 
— they generate shocks, hydrodynamic instabilities and turbulence in the hot intr- 
\Q acluster medium (ICM). They also cause "cold fronts", or contact discontinuities, 

CO around the cores of infalling subclusters or in the disturbed cool core of the main 

cluster. 1 In turn, these ICM disturbances distort and amplify magnetic fields (which 
are frozen into the ICM) and generate ultrarelativistic electrons that produce diffuse 
synchrotron emission of cluster-wide radio halosPEl Of these phenomena, with the 
• • current instruments, we can directly observe shock fronts and cold fronts as sharp 

. ^ brightness and temperature edges in high angular resolution X-ray images — Fig. [I] 

^ shows a poster child of the merger shocks, the "bullet cluster" IE 065 7-56 and 

Fig. [2] shows a cold front in the "cannonball galaxy" NGC 1404. 8 In the radio, the 
synchrotron emission is a product of density of ultrarelativistic electrons and B 2 , 
where the magnetic field B is deduced by indirect means. 9 Disentangling the rela- 
tivistic electron density and the magnetic field (e.g., by observing inverse Compton 
emission from the radio halos), as well as direct observation of turbulence in the 
ICM, are beyond the current observing capabilities. 

Figure [3] schematically shows the ICM density, pressure and specific entropy 
profiles in a sector in front of the "bullet" in the Bullet cluster. The sector crosses 
a cold front (the boundary of the bullet) and the shock that propagates ahead of 
it (Fig. [I]). While both these X-ray brightness edges exhibit similar density jumps, 
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(a) (b) 

Fig. 1. The Bullet cluster IE 0657-56. (a) Chandra X-ray image (0.8-4 keV band, 7' size), (b) 
gas temperature map (color; blue is T < 6 keV and yellow is T > 20 keV) and X-ray brightness 
(contours). A cool gas bullet is preceded by a hot shock front. 




sharp NW boundary of the galaxy is a cold front. 



the pressure across the cold front is almost continuous, as expected for a contact 
discontinuity between the dense, low-entropy remnant of an infalling subcluster and 
the surrounding shock- heated gas. The shock front exhibits a large pressure jump 
and a modest entropy increase, corresponding to its relatively low Mach number, 
M = 3. 

Cold fronts are more easily observable than shock fronts - they were first dis- 
covered with Chandra in a merging cluster A2142 10 and subsequently observed in 
most mergers. More surprisingly, they are also seen in the cool cores of at least half 
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Fig. 3. Schematic gas density, pressure and specific entropy profiles in a sector of the Bullet 
cluster crossing the bullet nose and the shock front. The density jump around r = 12" is a cold 
front (the front boundary of the gas bullet), and the jump at r = 50" is a shock front. 

of all relaxed clusters UHil Merger shock fronts quickly move away from the central, 
bright cluster regions into the faint outskirts, where they are difficult to observe in 
X-rays. To be discernible in an X-ray image, they also require a merger to occur 
almost exactly in the plane of the sky. For these reasons, only a handful is currently 
known, and they are the subject of this review. There are other types of shock 
fronts expected in clusters — far in the outskirts, a strong shock should separate 
the virialized cluster region and the outside infalling matter, but the X-ray bright- 
ness at those distances is far too low for these shocks to be observable at present. 
Explosions of the central AGNs also generate shocks in the cluster cores, usually 
with M « 1, but sometimes much stronger.^ These shocks are outside the scope of 
this review. 



2. Constraints on dark matter from the Bullet cluster 

We will start with some physical results already derived using the most promi- 
nent cluster shock front, the one in the Bullet cluster IE 0657-56. Fig. [I] shows an 
X-ray image and projected temperature map derived from a 500 ks Chandra ex- 
posurdUHlHI of IE 0657-56, and Fig. [Z] overlays a weak lensing total mass ma 
on the Chandra image. It is clear from these figures that IE 0657-56 is a merger 
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Fig. 4. X-ray image of IE 0657-56 (same as in Fig. [TJl) with contours showing total projected 
mass from weak lensing. The gas subcluster trails its mass peak because of ram pressure. This 
offset is direct proof of the existence of dark matterP^ 



of two subclusters flying apart just after a collision. The cooler gas bullet, with its 
characteristic shuttlecock shape, appears to be a remnant of the cool core of the 
smaller subcluster, pushed back from its host dark matter peak by ram pressure. 
The X-ray brightness and temperature profiles across the shock front give a Mach 
number (M = 3.0 ± 0.4) and a velocity of the shock front (4700 km s _1 )P The 
velocity of the subcluster is probably lower than that of the shock front, because 
the pre-shock gas flows in from the cluster outskirts toward the shock front under 
the gravitational pull of the subcluster. 17 The two subclusters should have passed 
through each other just 200 million years ago. 

It has been quickly realizecU^il^l that the offset between the peaks of the to- 
tal mass and the gas mass peaks seen in Fig. [4] offers the first direct and model- 
independent evidence for the existence of dark matter, as opposed to some forms 
of modified gravity, where the visible matter is all there is and the laws of gravity 
are incorrect on large scales 

18 19 

— an alternative possibility put forward to explain 
the longstanding problem of "missing mass" in galaxy clusters. 20 Indeed, the X-ray 
emitting gas is by far the dominant visible mass component in clusters, and yet, as 
the lensing mass map shows, the peaks of the total mass density are clearly located 
elsewhere. The measured total masses in those peaks are as expected from the ra- 
tio of gas mass, stellar mass in the galaxies, and total mass normally observed in 
clusters. What makes this cluster unique is the moment at which it is caught by the 
observer — when gas and dark matter have been temporarily spatially separated 
by a violent merger, revealing that these are indeed two different kinds of matter. 

The survival of two dark matter subclusters after a near-direct collision (in- 
dicated by the X-ray data) also provides an upper limit on self-interaction cross- 
section for the dark matter particles. The dark matter peaks coincide with peaks 
of the galaxy number density, and both dark matter and galaxy peaks are lo- 
cated ahead of the gas, as expected for collisionless dark matter and galaxies vs. 
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fluid-like ICM. The observed mass-to-light ratio within the gas-depleted subclus- 
ter peaks is found to be similar to that in other clusters within similar radii. This 
excludes the possibility that a significant fraction of dark matter particles has es- 
caped the subclusters as a result of particle collisions during the subclusters' recent 
passage through each other. Taking into account the merger velocity derived from 
the X-ray, this observation places a limit of a/m < 0.7 cm 2 /g on the dark matter 
self-interaction cross-section, excluding most of the astrophysically interesting 
range of the cross-sectionP^ 



3. Electron-proton equilibration timescale in the ICM 

The passage of a merger shock in a fully ionized intracluster plasma should heat the 
protons dissipatively, while electrons (for shocks with M <C (ra^/me) 1 / 2 « 43, which 
is always true for cluster mergers) are compressed adiabatically and subsequently 
heated by other mechanisms, such as Coulomb collisions with hotter protons. The 
linear sizes and temperatures of galaxy clusters are such that on a timescale of 
collisional electron-proton equilibration, r ep , the shock travels significant distances, 
creating observable regions where the electron temperature T e (the only one that 
we can currently measure in X-rays) is below the true thermodynamic tempera- 
ture. This is illustrated in Fig. |5j which shows simulated maps of the average and 
electron temperatures for a merging cluster with two shocks propagating outwards, 
assuming collisional r ep . The T e /T p nonequilibrium is expected in may types of as- 
trophysical objects, from solar wind to supernovae to WHIM filaments. However, 




(a) (b) 

Fig. 5. Simulations of a cluster merger,^ with colors showing the temperature and contours 
showing the mass. There are two shock fronts propagating ahead of two subclusters that have just 
crossed each other, (a) Average plasma temperature, (b) electron temperature assuming Coulomb 
electron-proton equilibration timescale. The linear scale is in Mpc and the temperature scale is in 
keV. The X-ray measured T e may be a significant underestimate of thermodynamic temperature 
at shocks. 
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Fig. 6. Predicted T e profiles for the shock in IE 0657-56, for the Coulomb electron-proton equili- 
bration timescale (blue band) and for instant equilibration (yellow band). Overlaid is the Chandra 
measurement (deprojected, with la error bars), which indicate a shorter equilibration timescale@ 

it is usually impossible to measure T e and T p independently and map them on 
linear scales on which their equilibrium is expected to be achieved. Shock fronts 
in clusters provide a unique opportunity to do so and thus constrain r ep . Because 
cluster shocks are relatively weak, the easily measurable plasma density jump at 
the shock is sufficiently far from its asymptotic value (factor of 4 for a 7 = 5/3), 
which allows one to derive the Mach number of the shock from the density jump us- 
ing the Rankine-Hugoniot jump conditions, and predict the equilibrium post-shock 
value of the plasma temperature independently of the measured jump in T e . The 
pre-shock T e gives the sound speed, which, together with the shock density jump, 
gives the speed of the post-shock gas flow relative to the shock front. This can be 
used to predict the post-shock T e profile for various values of r ep . For the shock in 
IE 0657-56 with M = 3, the expected rise in post-shock T e for the case of Coulomb 
collisions can be spatially resolved by Chandra, as shown in Fig. [6j The measured 
values of the post-shock T e exclude the Coulomb timescale at a 95% significance, 
favoring a much shorter r ep P 



4. Other clusters with known shock fronts 

Until recently, only two merger shock fronts were known and confirmed by the X-ray 
temperature measurement — the IE 0657-56 discussed above and A520 shown in 
Fig. [7^1. The latter has a shock with M ~ 2^ its recent long Chandra observation is 
currently being analyzed to validate the above constraint on r ep . A520 also exhibits 
a similar offset between gas and dark matter!^ 

Recently, several more shock fronts and shock candidates were found in X-rays. 
A unique case of two fronts located on the opposite sides of the cluster has been 
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(a) 

Fig. 7. 



Chandra image of A520; the shock front with M = 2 is to the SW of the bright central 
structure!^ (b) Chandra image of A2146, showing two M ~ 2 shocks about 300 kpc SE and NW 



of the center 



reported in A2146, 27 shown in Fig. [7p — a geometry expected for a symmetric 
merger with a small impact parameter. Both shocks have M ~ 2. Similar second 
shocks are not observed in either IE 0657-56 or A520, most likely because those 
shocks have already moved out of the bright central regions in these mergers of 
subclusters of very different masses. 

Now that we know how the merger shocks look like, it is possible to find them 
in archival X-ray data of lower quality. A possible shock front was detected in 
the ROSAT PSPC image of A75^l and recently confirmed using the Chandra 
temperature measurement across the shockpSlas shown in Fig. [8] This is a relatively 
weak shock with M = 1.6. 

As of this writing (August 2010), these are the only known merger shocks with 
the X-ray temperature profiles confirming that an edge-like feature in the X-ray 
brightness is indeed a shock front (and not a cold front, for example). There are 
several other likely shock fronts seen in the X-ray images but awaiting temperature 
confirmation. One is in Coma cluster, seen in a ROSAT PSPC mosaic (Fig. [9J. 
There is some indication of a temperature jump of the correct sign in the XMM- 
Newton temperature map, 31 though a measurement more matched to the shock 
position is needed to verify this. Coma is big enough for this putative shock front to 
be seen by the Planck observatory that maps the Sunyaev-Zeldovich signal, which 
is proportional to gas pressure and thus should readily reveal shocks. 

Another candidate front is found in the Chandra image of A2744, a spectacular 
merger at z = 0.3. An X-ray brightness profile in a sector containing the putative 
shock is shown in Fig. [10] The putative shock is a very low-contrast feature, likely 
corresponding to a low Mach number. We will see below that both Coma and A2744 
exhibit edge features in their radio halo maps that coincide with these putative 
shocks. 
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Fig. 8. (a) Chandra image of A754 with the sector showing a shock front with M ~ 1.6l 28 l 30 J (b) 
Chandra gas temperature profile showing a temperature jump across the front in a sector shown in 
panel (a) (radius is from the shock's center of curvature), confirming that this is indeed a shock. 30 




- 1 

******** ♦ 


1 1 1 - 




H 




: 




k 


1 


\ 

1 . tt : 



(b) 

Fig. 9. (a) ROSAT PSPC mosaic of Coma (1.6° x 1.6°), showing a brightness edge in the eastern 
sector (shown by dashes) that might be a shock front, (b) X-ray brightness profile in that sector, 
showing this feature at around r = 33 ; (red dash). 



It has long been suggested that peripheral radio relics are caused by electrons 
accelerated at shock fronts in the periphery of merging clusters. 32 However, most 
relics are located far from the X-ray bright cluster regions, so it is rarely possible to 
check for the presence of a shock front in an X-ray image, and even more difficult 
to get a confirming temperature measurement. 

A likely shock in A521 was discovered by analyzing the Chandra X-ray image at 
the position of a prominent radio relief see Fig. 11 There is a subtle X-ray bright- 
ness edge right where the "front" side of the relic is. The relic has a good radio 
spectrum, well represented by a power law. If one assumes that the radio-emitting 
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Fig. 10. (a) Chandra image of A2744 (10' x 10'), showing a subtle brightness edge in the SE 
sector (shown by dashes) that might be a shock front, (b) X-ray brightness profile in that sector, 
showing this feature at around r = 190" (red dash). 




(a) (b) 

Fig. 11. (a) Chandra image of A521, with GMRT 327 MHz synchrotron radio relic contours 
overlaid, 33 delineating a possible shock front, (b) Chandra X-ray brightness profile in the sector 
covered by the relic; radius is from the relic's center of curvature. Overlaid is a model shock profile 
with M = 2.3 predicted from the radio spectrum of the relief 3 - 3 -! which is consistent with the X-ray 
imaging data. Gas temperature data with sufficient quality is not yet available to confirm this 
shock. 



relativistic electrons come from Fermi acceleration at the shock, the spectrum im- 
plies M = 2.3. The density jump for such a shock is in good agreement with the 
X-ray brightness profile, as shown by red model in Fig.fTTjp. At lower radio frequen- 
cies, this cluster reveals a giant radio halo that starts at the relic and spans the 
entire cluster^ (Fig. 12 , discussed below). 



An X-ray brightness edge that looks like a shock was discovered in an XMM- 
Newton observation of the famous NW radio relic in A3667; 35 if this is indeed a 
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shock, the gas density jump indicates M ~ 2. A possible M ~ 2 shock front was 
discovered at the position of a radio relic in the cluster RXJ 1314-25P^The curious 
M shape of this front indicates a velocity gradient in the pre-shock gas — possibly 
an extreme case of the effect seen in the hydrodynamic simulations of the Bullet 
cluster P3 Table 1 gives a summary of the currently known merger shocks. 

Table 1. X-ray merger shock fronts and candidates (as of summer 2010) 



Cluster 



p jump T jump M Radio edge? X-ray refs. 



IE 0657-56 
A520 
A754 
A2146 N 
A2146 S 



yes 
yes 
yes 
yes 
yes 



yes 
yes 
yes 
yes 
yes 



3 
2 

1.6 

2 
2 



yes 

yes 

yes 
no data 
no data 



27 
27 



A521 

RXJ 1314-25 
A3667 
A2744 
Coma 



yes 
yes 
yes 
yes 
yes 



yes 
yes 
yes 
yes 
yes 



33 
36 
35 

this work 
this work 



5. Radio halos and X-ray shock fronts 

Merger shock fronts are expected and, as we saw above, observed to have relatively 
low Mach numbers. Indeed, a test particle falling into the gravitational potential of a 
typical cluster would acquire M ~ 3 — 4. An infalling subcluster would generate mul- 
tiple shocks in front of it, preheating the ambient gas and further reducing its Mach 
number by the time it arrives into the cluster's X-ray bright central region, where we 
observe it. Such weak shocks are believed to be inefficient accelerators of electrons 
from the thermal pool to Lorentz factors ~ 10 4 required to produce relics and radio 
halos E2E3 Yet, it has been noticed that synchrotron radio halos in IE 0657-5(P and 
A52CPS1 display brightness edges coincident with the X-ray shocks (Fig. 12), indi- 



cating that these shocks must have something to do with producing radio emission 
at least at those locations — either accelerating electrons, or compressing or re- 
accelerating a pre-existing relativistic electron population, and perhaps strengthen- 
ing the magnetic fields beyond simple compression.^ Re- acceleration of pre-existing 
relativistic electrons is an attractive theoretical possibility! 38 ! 39 ! 

All the clusters with newly discovered X-ray shocks and shock candidates exhibit 
radio halos or relics, except A2146 that doesn't have sensitive radio data yet. It is 
interesting to overlay the radio images on the X-ray shocks, which we do in Fig. [l2j 
In addition to IE 0657-56 and A520 discussed above, it shows A521, which exhibits a 
relic coincident with the probable shock, plus a halo extending from that relic across 
the whole clusteJUMI — a radio morphology very similar to A520. A754, with the 
weakest shock (M = 1.6), exhibits a steep-spectrum relic at 74 MHz SSI coincident 
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(A2744) (Coma) 

Fig. 12. Contours o f rad io halos overlaid on X-ray images of clusters with shock fronts 
(IE 0657-56[M3 A520(2£|43J or front candidates (A521, 34 A2744, 44 Coma 45 ). In all cases, 

there is coincidence of a shock front seen in the X-ray with an edge in the radio halo. 
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with the shock front, and a radio halo extending from the relic westward across the 
clusterESES Finally, the two shock candidates in A2744 and Coma, which we report 
in this paper, both have edge-like features seen in recent high-quality radio images 
of their halosJHES coincident with putative X-ray shock fronts. In Coma, the X-ray 
edge spans the southern half of the radio edge (Fig. [l2|. The cluster RX J 1314-25 
(not shown) also has a relic at the shock front and and a small radio halo extending 
from itP3 

Thus, spatial coincidence of the merger shocks with edges of radio halos — or 
with radio relics that delineate an edge of the radio halo — is quite ubiquitous. 
It supports the scheme that we proposed for A520 and A52lj^ in which a weak 
merger shock re-accelerates pre-existing relativistic electrons (which may remain 
from earlier mergers, or generated by collisions of the long-lived cosmic ray protons 
with thermal protons^). As these electrons move downstream from the shock, they 
rapidly cool, creating a narrow arc-like "relic" , but at a certain distance are picked 
up and re- accelerated again by turbulence generated by the merger behind the 
shock, which produces the cluster-wide halo. In this picture, the edge of the radio 
halo and the bulk of the halo are distinct phenomena, though both caused by the 
same merger. Depending on the Mach number of the shock, the magnetic field 
behind the shock (that determines the rate of synchrotron cooling of the electrons 
and thus the width of the "relic" in the direction across the shock front), and the 
power spectrum of turbulence in the body of the ICM behind the shock, the relic 
may appear at some frequencies as a distinct radio source, while at other frequencies, 
merge with the halo. In A521, the relic is dominant at v > 1 GHz, while in A754, 
the relic dominates at 74 MHz but merges with the halo at higher radio frequencies. 
The prediction of this picture is that the radio spectrum at the halo "edge" should 
be a power law determined by the shock's Mach number (Fermi acceleration from 
thermal pool and re-acceleration from plausible initial electron distributions result 
in the same spectrum). This is again consistent with observations of A521^ and 
A754pj rac Q S p ec trum of the rest of the halo (more exactly, the frequency 
of its exponential cutoff) is determined by the velocity of the turbulence and the 
strength of the magnetic field .^Depending on that spectrum, it may or may not be 
possible to see the steepening of the radio spectrum of the front edge due to aging 
of electrons as one moves from the shock front inward (as was possible for A521^). 

We note that some more irregularly shaped cluster radio relics, such as a bright 



relic in A2744 (seen in the left edge of the image in Fig. 12), A2256 or A3667, may 
have different origin, perhaps involving a shock passage across a distinct region of 
fossil radio plasmaP^ 

It is clear that detailed, spatially resolved studies of cluster merger shocks in the 
X-ray and at multiple radio frequencies is a promising way to study the cosmic ray 
acceleration mechanisms in astrophysical plasmas. We may be starting to collect a 
sufficiently large sample for such studies. 
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